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Young stellar population and ongoing star formation in the 
H II complex Sh2-252 
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ABSTRACT 

In this paper, an extensive survey of the star forming complex Sh2-252 has been un- 
dertaken with an aim to explore its hidden young stellar population as well as to 
understand the structure and star formation history for the first time. This complex 
is composed of five prominent embedded clusters associated with the sub-regions A, 
C, E, NGC 2175s and Teu 136. We used 2MASS-NIR and Spitzer-IRAC, MIPS pho- 
tometry to identify and classify the young stellar objects (YSOs) by their infra-red 
(IR) excess emission. Using the IR colour-colour criteria, we identified 577 YSOs, of 
which, 163 are Class I, 400 are Class II and 14 are transition disk YSOs, suggesting a 
moderately rich number of YSOs in this complex. Spatial distribution of the candidate 
YSOs shows that they are mostly clustered around the sub-regions in the western half 
of the complex, suggesting enhanced star formation activity towards its west. Using 
the spectral energy distribution (SED) and optical colour-magnitude diagram (CMD) 
based age analyses, we derived probable evolutionary status of the sub-regions of Sh2- 
252. Our analysis shows that the region A is the youngest (~ 0.5 Myr), the regions B, 
C and E are of similar evolutionary stage (~ 1-2 Myr) and the clusters NGC 2175s and 
Teu 136 arc slightly evolved (~ 2-3 Myr). Morphology of the region in the 1.1 mm map 
shows a semi-circular shaped molecular shell composed of several clumps and YSOs 
bordering the western ionization front (IF) of Sh2-252. Our analyses suggest that next 
generation star formation is currently under way along this border and that possibly 
fragmentation of the matter collected during the expansion of the H n region as one of 
the major processes responsible for such stars. We observed the densest concentration 
of YSOs (mostly Class I, ~ 0.5 Myr) at the western outskirts of the complex, within 
a molecular clump associated with water and methanol masers and we suggest that 
it is indeed a site of cluster formation at a very early evolutionary stage, sandwiched 
between the two relatively evolved CH II regions A and B. 

Key words: stars: formation — stars: pre— main— sequence — infrared: ISM — 
H II regions - ISM: individual objects: Sh2-252 



1 INTRODUCTION various star forming regions at different environments are 

, , „ , . , very essential. Census of young stellar objects (YSOs) born 

In order to understand the interstellar processes which gov- . , , ., . , , . , ,. , ., .: , .. 

ma cloud, their mass, age and spatial distribution, together 
ern star formation, studies of the large scale properties ol ... ,,, , . , , ,,, , , , 

° with mass ol the star forming cloud are some of the key data 

used for the assessment of the time scale and star formation 
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ber stars of a given young association is quite difficult, ft 
is well known that YSOs show excess emission in infra-red 
(IR) above main sequence (MS) photosphere due to ther- 
mal emission from their circumstellar material. Thus, YSOs 
can be identified by looking for IR excess emission. YSOs 
are oft en categorized into Cla ss 0, I, II or III evolutionary 
stages (|Lada fc Wilking||l984 ). Class objects are deeply 
embedded protostars that are still experiencing cloud col- 
lapse. They are extremely faint at wavelengths shorter than 
10 /J,m and have a significant sub-millimeter luminosity. A 
Class I YSO is also an object whose emission is dominated 
by a dense infalling spherical envelope, however these ob- 
jects are bright in the IR. A Class II YSO is characterized 
by the presence of an optically thick, primordial circum- 
stellar disk, which dominates the star's emission. When the 
circumstellar disk material becomes optically thin, the star 
is classified as a Class III star. Hence, IR survey of a star 
forming region can be used as a powerful tool to distinguish 
the stars with IR excesses from stars without such excesses. 
The unprecedented sensitivity and mapping capabilities of 
the Spitzer Space Telescope provides an excellent platform 
to survey the star forming regions in the mid-IR and to 
identify the IR excess sources in the region. Furthermore, 
Spitzer data can also be used to classify YSOs at different 
evolutionary stages, which helps us pro bing the recent sta r 
formation activity of a given region fe.g.. lSmith et al.ll2010l ). 

Massive stars have a profound effect on their natal en- 
vironment creating wind-blown shells, cavities and HII re- 
gions. The immense amount of energy released through their 
stellar winds and radiation disperse and destroy the remain- 
ing molecular gas and likely inhibit further star formation. 
However, it has also been argued that in some circumstances, 
the energy input by these massive stars can promote and in- 
duce subsequent low mass star formation in the surrounding 
molecular gas before it disperses (e.g., Koenig et al. 2012). 
Therefore, regions that contain massive stars provide excel- 
lent laboratories for studying both high and low mass star 
formation. Identification and characterization of the YSOs 
in star forming complexes hosting massive stars are essen- 
tial steps to examine the physical processes that govern new 
generation star formation in such complexes. 

Here we present the results of our multi- wavelength 
analy ses of the star forming complex Sh2-252 l|Sharplessl 
1959; Q2ooo = 06 h 09 m 39 s ; 8 20l 



+20°29'12"; 1=190°. 04; 



b=+0°.48). Sh2-252 is an evolved H II region powered by the 
central massive ionizing source of 06.5V type, HD 42088, 
which is a member of the Gemini OBI association. The 
DSS2-R band image of the region around Sh2-252 for an 
area of ~ 1.3 x 1.3 deg 2 is shown in Fig. Q] Radio obser- 
vations at 5 GHz detected six extended sources (Sh2-252 
A to F) at different spatial locations towards the region, of 
which, four of them (A, B, C and E) are classified as com- 
pact H ii (CH ii ) regions JFelli et al.ll 19771 . lLada fc Wooden! 
1 19791 ) . The stellar cont ents of this region were not so well 
explored until recently. Ijose et al.l (|2012h (hereafter Paper 
1) unraveled its stellar contents using deep optical, shallow 
near-infrared (NIR; 2MASS), Ha survey along with spectro- 
scopic observations. They estimated the distance (2.4 kpc), 
reddening and identified 12 OB stars of spectral type ear- 
lier than B8, 61 Ha emission line sources within the nebula 
and also 5 prominent embedded clusters associated with the 
sub-regions A, C, E, NGC 2175s and Teu 136. The age and 



age spread, initial mass function (IMF), K-band luminosity 
function (KLF) etc. of the sub-regions have also been dis- 
cussed in Paper 1. The main sub-regions of Sh2-252 such 
as A, B, C, E, F, NGC 2175s and Teu 136 are marked 
in Fig. [T] Sh2-252 is bounded by an ionization front (IF) 
on t he west side, wh ile it is density bounded on the east 
side (|Felli et al.lll977fl. The detailed 12 CO and 13 CO maps 
of the region by lLada fc Wooden! (|1979T ) showed that the 
cloud complex is separated into two fragments by a long 
rift of little CO emission and named them as western cloud 
fragment (WCF) and eastern cloud fragment (ECF), respec- 
tively, and are marked in Fig. [T] Most of the mass of WCF 
is concentrated in a narrow ridge which borders the IF of 
Sh2-2 52 at its west (JLada fc Wooden! 1 19791 ; iKoempe et"al] 
Il989f ). The similarity in the kinematics of these two cloud 
fragments shows that they belong to the same cloud complex 
and ECF probably lies at the back edge of the H II region, 
though it is projected very close to the i onizing source 
(|Lada fc Woodenlil979l ; iFountain et alj|l983l ). The most in- 
tense CO peak of the WCF is located very close to region A, 
with water and methanol m a ser emissions (|Lada fc Wooden! 
1 19791 ; IKoempe et~aH Il989l ; ISzvmczak et al.1 l200d ) within 
its proximity, sug gesting recent star formation activity to- 
wards this region (|Lada fc Woodenlll979l ; lLada et aUll98ll ; 
IKoempe et al.l ll989). An extensive overview of the previous 
studies of Sh2-252 can be found in Paper 1. 

This work is the continuation of Paper 1, where we 
have focused mainly on the characteristics of the YSOs in 
the Sh2-252 complex for the first time and to constrain the 
star formation history of the complex. We used the optical, 
2MASS-JHK, Spitzer- IRAC and MIPS data sets to iden- 
tify and classify the YSOs of Sh2-252 region. An attempt 
has also been made to characterize the physical properties 
of the YSOs based on the spectral energy distribution (SED) 
analyses. Section 2 describes various data sets used for the 
present study. Analyses and results including the identifica- 
tion and classification of YSOs, their various properties such 
as ages, masses, SEDs, spatial distribution etc. have been 
discussed in Section 3. Section 4 deals with discussion on 
the morphology and star formation scenario of the complex 
and finally summary of the paper is presented in Section 5. 



2 DATA SETS 

2.1 Spitzer-IRAC observations 

Two sets of IRAC observations in 3.6, 4.5, 5.8 and 8.0 
/un bands (channels 1, 2, 3 and 4) were available in the 
Spitzer space observatory archive program (Program IDs: 
201; 20506, Pis; G. Fazio; J. Hester) which were taken to- 
wards WCF and ECF on 2004 October 28 and 2006 March 
26, respectively. The IRAC data were taken in the High 
Dynamic Range mode with three dithers per map position 
and with integration time of 0.4 sec and 10.4 sec per dither. 
We downloaded the basic calibrated data images (versions 
S18.7.0 and S18.5.0) from the Spitzer archivqj and the raw 
data were processed and calibrated with the IRAC pipeline. 
The final mosaics were created using the MOPEX pipeline 
(version 18.0.1) with an image scale of l".2 per pixel. The 
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Figure 1. The DSS2-R band image of the r egion around Sh2 -252 for an area of ~ 1.3 X 1.3 deg 2 . The locations of the thermal 
radio sources A, B, C, E and F identified bv iFelli et alj q 1977T 1 ■ and the small clusters NGC 2175s and Teu 136 are marked in the 
figure. The ionizing source (HD 42088) o f the H II region i s mark ed by using a star symbol. The contours are the 12 CO Raylcigh-Jeans 
brightness temperature map taken from lLada &: Wooden! IJ1979T ). The contours are shown with an interval of 3 K, where the lowest 
contour corresponds to the temperature of 6 K. The region shown using the dashed box is covered by the Spitzer-IRAC observations 
and the regions shown by solid boxes represent the area covered for the Spiizer-MIPS-24 jim observations. The details of the WCF and 
ECF are given in the text. 



total area covered by these two observations is marked in 
Fig. Q] by using a dashed box. 

We performed the point response function (PRF) fitting 
method in multiframe mode using the tool APEX developed 
by Spitzer Science Center on all the Spitzer-IRAC images to 
extract the magnitudes. The standard PRFmap tablqj pro- 
vided on the Spitzer website were used to fit variable PRFs 
across the image. Point sources with peak values more than 
5cr above the background were considered candidate detec- 
tion. Many sources detected in the nebulosity were appeared 
to be spurious. The spurious sources were identified visually 
and deleted from the automated detection list. Also, some of 
the sources were not detected by APEX automatically. We 
identified these sources visually and added the coordinates in 
the user-list mode in APEX to extract the magnitudes and 



http://ssc.spitzer.caltech.edu/irac/calibrationfilcs/psfprf/prfmap.tbl 



thus we made sure that photometry of every genuine source 
is derived. We have adopted the zero-points for conversion 
between flux densities and magnitudes to be 280.9, 179.7, 
115.0 and 64.1 Jy in the 3.6, 4.5, 5.8 and 8.0 fj,m bands, re- 
spectively, following the IRAC Data Handbook. There were 
signs of saturation of bright sources in the long integrated 
images. The photometric data derived from the short im- 
ages were replaced for the saturated sources. The two sets 
of IRAC observations had substantial overlapping between 
the regions and we used the average magnitudes of those 
common sources detected in both the observations. Chan- 
nels 3 and 4 detected far fewer sources than channels 1 and 
2. This deficit exists because channels 1 and 2 are more sen- 
sitive than channels 3 and 4 and are also less affected by the 
bright diffuse emission that dominates the channel 3 and 4 
images. The underlying typical stellar photosphere is also 
intrinsically fainter at channels 3 and 4 than at channels 1 
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and 2, which further hinders the ability to detect sources in 
the longer wavelength IRAC filters. We finally obtained ex- 
traction lists of 24877, 20591, 5083 and 3263 sources in the 
IRAC 3.6, 4.5, 5.8, and 8.0 /jm bands, respectively, within 
the area shown in Fig. [l] The IRAC data of the four band 
passes were merged, by matching the coordinates using a ra- 
dial matching tolerance of 1".2. To ensure good photometric 
accuracy, we restricted our catalog with sources having un- 
certainty ^ 0.2 mag in all the bands. Thus our final IRAC 
catalog contains photometry of 31251 sources which were 
detected in one or more IRAC bands. Out of these, only 
1988 sources were detected in all the IRAC bands. 



2.2 Spitzer-MIPS 24 (im observations 

The areas towards WCF and ECF were observed at 24 
/xm using the Multi-band Imaging Photometer for Spitzer 
(MIPS) on 2006 October 13 and 2008 October 25 (Pro- 
gram IDs: 20726; 50758, Pis: J. Hester; M. M. Drosback). 
We downloaded the BCD images (version S18.12.0) from 
the Spitzer archive and the final mosaics were created using 
the MOPEX pipeline (version 18.0.1) with an image scale of 
2". 45 per pixel. The total area covered by these observations 
is marked in Fig. [l] The sources were identified visually and 
their coordinates were given in the user-list mode in APEX. 
To extract the flux, we performed the PRF fitting method 
in the single frame mode. The zero-point value of 7.14 Jy 
from the MIPS Data Handbook has been used to convert 
flux densities to magnitudes. The final catalog contains the 
24 fj,m photometry of 646 sources having uncertainty ^0.2 
mag. We cross-correlated our MIPS-24 /im catalog with the 
source catalog of IRAC, using a 2". 5 matching radius. We 
obtained 24 jj,m counterparts for 428 sources with IRAC 
data in all the bands. The 24 /im photometry was mainly 
used to classify the YSOs and to analyze their SEDs (see 
Section |3"XTJ) . 



2.3 Near- infrared data from 2MASS 

NIR JHK data for point sources within a radius of 30 ar- 
cmin around HD 42088 hav e been obtained from 2MASS 
point source catalog (PSC) (|Cutri et al.ll2003h . To improve 
the photometric accuracy, we used the data with photomet- 
ric quality flag (ph_qual = AAA) which gives a S/N ^ 10 
and a photometric uncertainty < 0.1 mag. This selection 
criterion ensures best quality detection in terms of photom- 
etry and astrometry as given on the 2MASS websitqj. We 
used 2MASS data along with Spitzer-IRAC data to identify 
and classify the candidate YSOs within Sh2-252 (see Section 
13.3.21) . The NIR and IRAC sources were bandmerged within 
l".2 radial matching tolerance. If more than one 2MASS 
sources satisfied this requirement, then the nearest one was 
selected. Thus we obtained a total of 9097 2MASS sources 
having IRAC counterpart at least in one band. 



3 RESULTS 

3.1 General morphology of the region 

A colour composite image using the 3.6, 4.5 and 8.0 /zm 
bands of Spitzer-IRAC is shown in Fig. [2] The important 
sub-regions and the ionizing source (HD 42088) of Sh2-252 
are marked in the figure. Different IRAC bands have contri- 
bution from various extended emission features. For exam- 
ple, 3.6 fim band is dominated by diffuse emission from the 
3.3 /im PAH feature. The 4.5 /im band includes the Bra line 
plus CO band head emission and Hi lines, and 5.8 and 8.0 
/xm bands are dominated by strong PAH features, although 
a few sources show extended thermal continu um emissions 
from warm circumstellar dust in these bands (jSmith et al.l 
|2010j). PAHs are believed to be destroyed in the ionized gas 
and thought to be excited in the photon dominated region 
(PDR) that is situated at the interface of H 11 region and 
molecular cloud by the ab sorption of far-UV p hotons leak- 
ing from the H n region (jPomares et al J 12009 ). Therefore 
PAHs are good tracers of the warm PDR that surrounds the 
H 11 region. The most striking feature in Fig. [2] is the ma- 
genta coloured emission that appears almost semi-circular 
in shape towards the west of Sh2-252. The PAHs within the 
PDR re-emit their energy in the 3.6, 5.8 and 8.0 (im bands, 
and this can be seen as magenta coloured emission in the 
colour composite image given in Fig. [2] Altogether, this im- 
age shows a bubble-like emission feature around the ionizing 
source HD 42088, likely to be created due to the expansion 
of the H 11 region to its surroundings. The bright-rim feature 
discussed in Section[T](i.e., region F; see also Fig. 2 of Paper 
1 for a close-up view) coincides well with the PAH emission 
feature seen at the southwest of Fig. [2] with a sharp edge 
pointing towards HD 42088. This sharp edge indicates the 
IF, that coincides with the molecular ridge seen in the 12 CO 
map shown in Fig. [l] as well as the extinction map shown 
in Fig. [3] (see Section [33J- This correlation indicates the in- 
terface of the ionized and molecular gas. The correlation of 
the half-ring structure of PAH emission towards the west of 
Sh2-252 with the 12 CO emission and extinction maps indi- 
cates that the ionized gas is separated from the molecular 
cloud by the neutral gas with PAH emission. The absence of 
8 fim emission in the interior of the H 11 region can be inter- 
preted as the destruction of PAH molecules by the intense 
UV radiation of the ionizing source. There is also less intense 
diffuse emission towards the east of Sh2-252. However, the 
ionized gas is bounded more extensively in the western di- 
rection than in the eastern direction. It is to be noted that 
the location of NGC 2175s is devoid of strong PAH emission 
features. Interestingly, another bright red emission knot is 
seen further towards east of NGC 2175s with a bright red 
point source inside it. A small clustering is evident further 
east to this point source which coincides with the sub-region 
Teu 136 discussed in Paper 1. There are not many studies on 
this group of stars in the literature. However, based on our 
spectral analysis in Paper 1, we concluded that this group 
is a part of the main Sh2-252 complex itself. 



http://www.ipac.caltcch.edu/2mass/releases/allsky/doc/ 




Figure 2. IRAC three colour image of Sh2-252 (blue: 3.6 (an; green: 4.5 fim; red: 8.0 /im; see the electronic version for the colour image). North is up and east is to the left. The 
important sub-regions are marked in the figure and the white star symbol represents the location of HD 42088, the ionizing source of Sh2-252. 
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3.2 A'-band extinction map from H and K data 

The indication of variable extinction in the Sh2-252 region 
has been discussed in Paper 1 based on the optical photo- 
metric analyses of the bright stars as well as from the distri- 
bution of candidate YSOs on the (J — H)/(H — K) colour- 
colour (C-C) diagram. In order to quantify the amount of 
extinction within each sub-region of Sh2-252 and to charac- 
terize the structure of molecular clouds, we derived the Ak 
extinction map using the (H — K) colours of the stars in 
the 2MASS catalog. To improve the quality of the map, the 
candidate YSOs and probable contamination sources (see 
Section I3.3p were excluded from the list. We used the grid 
method to determine the mean value of Ak- Briefly, we di- 
vided the Sh2-252 region into many grids of size 30" x 30", 
and the median value of (H — K) colours of all the sources 
within each grid was measured. The sources deviating above 
3cr were excluded to calculate the fi nal median colour of e ach 
grid. We used the reddening law bv lFlahertv et al.l (|2007f ) to 
convert (H — K) colour in to Ak, using the relation Ak — 
1.82 x {H - K) obs -(H- K) int , where, (H - K) in t = 0.2 
is assumed as an average intrinsic colour of stars in young 
clusters (see lAllen et al)l2008l ; rGutermuth et alj2009h . Since 
this complex is situated at a distance of 2.4 kpc (Paper 1), 
the contribution from the foreground field stars would have 
affected the present extinction measurements. An average 
foreground extinction of Av= 1.5 mag has been estimated 
towards this region in Paper 1. Therefore, to eliminate the 
foreground contribution in the extinction measurement, we 
used only those stars with Ak > 0.15 mag to generate the 
extinction map. There were an average of five sources in 
each grid to calculate the median colour. The final extinc- 
tion map, having an angular resolution of 30", which is sen- 
sitive down to Av = 16 mag is shown in Fig. [3] However, 
it has to be kept in mind that this map is limited by the 
sensitivity of the 2MASS survey. As a result, the Ak value 
derived may be underestimated due to the small number of 
background stars detected in the heavily extincted area of 
the complex. 

It is interesting to note that the extinction map derived 
from the (H — K) colours shown in Fig. [3] resembles the 
general distribution of the molecular cloud as outlined by 
the 12 CO emission map shown in Fig. [T] It traces the pres- 
ence of molecular cloud towards the sub-regions A, B, C and 
nicely follows the bright rim feature visible at the west of 
Sh2-252 (i.e., region F, see Fig. 2 of Paper 1) similar to the 
ridge of molecula r cloud seen in the 12 CO emission map by 
lLada fc Wooden! (|1979r ). Similarly, the separation between 
the two cloud fragments WCF and ECF (see Section [1} by 
the long rift of less extinction is also clearly seen in Fig. [3] 
The western cloud portion is more extincted as compared 
to the eastern one. However, isolated extinction complexes 
are seen in the eastern side as well. It has to be noted that 
the locations of the clusters NGC 2175s and Teu 136 are de- 
void of heavy reddening. The same trend has already been 
noticed in Paper 1. The average extinction Ay of the four 
optically bright members of the cluster NGC 2175s (see Ta- 
ble 3 of Paper 1) is ~ 2.2 mag, which is significantly lower 
than the extinction towards the western part of the com- 
plex. This could be explained by the fact that the strong 
stellar winds from the four early type sources in NGC 2175s 
would have removed the dust and gas from its neighborhood 



or the region could be more evolved than the other parts of 
the complex. 



3.3 Identification of YSOs based on near- and 
mid-infrared colours 

In the ensuing sections, we identify and classify the YSOs 
within Sh2-252 mainly into Class I and Class II categories 
using their IR colours. We do not classify diskless Class III 
YSOs because they are indistinguishable from field stars 
based on their IR colours and hence we cannot reliably iden- 
tify them using the current data set. The main limitation 
of the identification of YSOs based on their IR colours is 
the contamination from various non-stellar sources in IRAC 
detection such as, polycyclic aromatic hydrocarbon (PAH) 
emitting galaxies, broad-line active galactic nuclei (AGNs), 
unresolved knots of shock emission, PAH-emission contam- 
inated apertures etc., which have colours similar to that of 
YS Os and may lead to wr ong id entification o f YSO s. 

iMegeath etail (|2004h and lAllen et al.l (|2004h intro- 
duced the classification criteria for YSOs using the ([3.6] — 
[4.5], [5.8] - [8.0]) C-C diagram. These colour crite- 
ria si nce then have bee n mod i fied by several authors 
(e.g lHarvev et al.l I2006L 120071 . IGutermuth et all 120081 . 
12003 . I Evans et al.l 12009 and the references therein). We 



adopted the emp i rical three-phase scheme introduced by 
IGutermuth et al.l (|2009h to identify the YSOs with IR ex- 
cess in this complex. Briefly, in phase 1, after excluding the 
non-stellar contamination (see Section l3.3.ip . all four IRAC 
bands are used to separate the sources with IR excesses due 
to disks and envelopes (Class II and Class I, respectively) 
from pure photospheres (Class Ill/field stars) based on [3.6] 
- [5.8] and [4.5] - [8.0] colours. Phase 2 is applied to those 
sources which lack IRAC detection at [5.8] or [8.0] /im, but 
which have high-quality 2MASS detection (a < 0.1 mag) in 
at least H and K bands. Sources are then classified as Class 
II or Class I objects based on their de-reddened (([K] — [3.6] )o 
and ([3.6] - [4.5])o) colours (|Allen et alj|2008h . In addition, 
there are a few objects with NIR+IRAC colours of Class 
Ills, but which exhibit a 24 [im excess. These objects are in- 
cluded in a new class of YSOs i.e., 'transition disk' sources. 
Phase 3 re-examines all sources including contaminants and 
identifies Class I objects and transition disk objects on the 
basis of their MIPS 24 (im photometry. 



3.3.1 Selection of YSOs using IRAC data 

In order to identify the YSOs by their IRAC colours, we 
used only those sources detected in all the IRAC bands 
with photometric uncertainties less than 0.2 magnitude. 
Various non-stellar contaminating sources mentioned above 
are eli minated by using severa l criteria described in de- 
tail by IGutermuth et al.1 (|2009h . Out of the 1988 sources 
detected in all the IRAC channels, 444 sources are found 
to be contaminants. The remaining sample of 1544 sources 
was used to identify the Class I and Class II YSOs associ- 
ate d with the region using the IRAC colour criteria given 
by IGutermuth et al.l (|2009i l. The above criteria yields 340 
YSOs with IR excess emission, with 86 and 254 sources 
having colours consistent with Class I and Class II types, re- 
spectively. The rest could be Class Ill/field stars, which we 
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Figure 3. A"-band extinction map generated from the (H — K) colours of 2MASS data. The contours begin at Ak = 0.55 mag and 
increase by 0.15 mag up to 1.6 mag. The locations of the important sub- regions of Sh2-252 are marked in the figure. The abscissa and 
the ordinate are for the J2000 epoch. 



are unable to classify reliably based on their IRAC colours. 
However, it is to be noted that in a few cases, highly red- 
dened Class II sources could m imic the colours of Class I 
source (|Gutermuth et al J 12003 ). In Fig. [4] (left panel), we 
have shown the [3.6] - [4.5] vs [5.8] - [8.0] C-C diagram for 
all the uncontaminated IRAC sources, where the Class I and 
Class II sources are shown as red squares and green trian- 
gles, respectively. Here we would like to point out that the 
present YSO selection based on the IRAC- four colour de- 
tection is not complete, as many sources falling in the bright 
nebulous regions were not detected in the 5.8 and 8.0 lira 
bands. 



3.3.2 Additional YSOs from H, K, 3.6, 4.5 and 24 /im data 

The detection of Class I and Class II YSOs in all the four 
IRAC bands is limited mainly by the lower sensitivity and 
enhanced nebulosity of the [5.6] and [8.0] fira bands. In or- 
der to account for the YSOs missing in these bands, we used 
H, K, [3.6] and [4.5] /im data sets and the addi t ional YSOs 
in the complex are identified. ICutermuth et al.l 1)2003 ) used 
the dereddened colour criteria to identify such YSOs. We 
used the extinction map (see Fig. [3] Section I3.2f) to dered- 
den the data. To minimize the inclusion of extra-galactic 



contaminants, we applied a simple brightness limit in the 
dereddened [3.6] /ira photometry i.e., all the Class I YSOs 
must have [3.6]o ^ 15 mag and all the Class I I YSOs must 
have [3.6]o < 14.5 mag (|Gutermuth et alj|200a). Thus, after 
remov ing the contaminants and using the lCutermuth et al.l 
l|2009f ) colour criteria, we identified a total of 367 candidate 
YSOs (67 Class I; 300 Class II) with IR excess emission, of 
which, 211 (52 Class I and 159 Class II) are found to be com- 
mon with the YSOs identified using the IRAC four-colour 
criteria (see Section [3]3T]) . Finally, we obtained 15 and 141 
additional YSOs having colours consistent with Class I and 
Class II, respectively using the H, K, [3.6] and [4.5] colour 
combination. In Fig. [4] (middle panel), we have shown the 
{[K] - [3.6]) vs ([3.6] - [4.5])o C-C diagram for all the un- 
contaminated sources, where the Class I and Class II sources 
are shown in red squares and green triangles, respectively. 

In order to identify the transition disk sources, we re- 
examined the entire catalog of sources having 24 ^m coun- 
terparts and with uncertainty less than 0.2 mag. Those 
sources that were considered photospheric in the previous 
classification but have significant excess emission at 24 fira 
(i.e., [5.8]-[24] > 2.5 mag or [4.5]-[24] > 2.5 mag) are typ- 
ically known as transition disks. They are the Class II 
sources with significant dust clearing within their inner disks 
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Figure 4. left: IRAC C-C diagram of all the uncontaminatcd source s identified within the re gion shown by a dashed box in Fig. [T] The 
YSOs classified as Class I and Class II, based on the colour criteria bv lGutermuth et alj 1 120091) , are shown by using red squares and green 
triangles, respectively. The reddening vector for Ak = 5 mag is plotted by using the reddening law from iFlahertv et al.l 1120071 ). middle: 
The dereddened ([K] — [3.6])o vs ([3.6] — [4-5])o C-C diagram for all the IRAC sources having 2MASS counterparts after removing the 
contaminants. The sources in red and green are the candidate Class I and Class II YSOs, respectively, right: [3.6]-[4.5] vs [8.0]-[24] C-C 
diagram for all the IRAC sources having 24 /im photometry after removing the contaminants. The Class I, Class II and transition disk 
YSOs are shown by using red squares, green triangles and magenta circles, respectively. 



(|Gutermuth et al.M2009f ) and we identified 14 such sources 
within the complex. Also, any source which lacks detection 
in some IRAC bands yet has bright 24 ^m photometry are 
likely to be deeply embedded protostars. We identified 67 
such sources and those are included in the Class I category 
in the following analysis. Five sources which lacked strong 
24 /im detection but had been classified as Class I YSOs in 
the previous step (Section I3.3.ip are reclassified as heavily 
reddened Class II YSOs. 

Thus after reclassifying the sources using the 24 ^m 
photometry, we have finally identified 577 YSOs with IR 
excess emission in the Sh2-252 complex. Of these, 163 are 
consistent with Class I, 400 are consistent with Class II and 
14 are consistent with transition disk YSOs. These objects 
will be used in the following analysis. Thus the Sh2-252 com- 
plex seems to have moderately rich number of YSOs. In the 
right panel of Fig.g]we have shown the [3.6]-[4.5] vs [8.0]-[24] 
C-C diagram of all the uncontaminated IRAC sources hav- 
ing counterparts in 24 /jm along with the candidate Class I, 
Class II and transition disk sources marked in red squares, 
green triangles and magenta circles, respectively. 



3.4 Completeness of the census of YSOs 

To evaluate the completeness of the census of IRAC detec- 
tion quantitatively, we plot histograms of the IRAC sources 
having error < 0.2 mag for all the channels and they are 



shown by using solid lines in Fig. [5] In general, the data can 
be considered complete till the linear distribution in the his- 
tograms which are found to be 16.5, 16.0, 14.0 and 13.0 mag 
for 3.6, 4.5, 5.8 and 8.0 ^m bands, respectively. However, the 
completeness of the YSO census is limited by many factors. 
Bright extended, variable nebulosity in the IRAC bands can 
be found across the region, which significantly limits the 
point source detection in these areas. The YSO identifica- 
tion from the 2MASS- IRAC colours is limited by the modest 
sensitivity of the 2MASS survey, while that from the IRAC- 
MIPS colours suffer from the significant saturation in the 
IRAC 8 /im and MIPS 24 /jm images caused by the central 
luminous sources as well as bright nebulosity. Similarly the 
variable reddening (see Fig. [3| and stellar crowding charac- 
teristics across the extend of the region (see Paper 1) could 
also affect the local completeness limit. Completeness lim- 
its for the Spitzer photometry were also assessed by using a 
method of inserting artificial stars into the mosaics and then 
employing our d etection algorithms to identify them (e.g., 
I Jose et alj 120081 ). The sub-regions A, B, C, E and F of the 
Sh2-252 complex are found to be severely affected by the 
variable reddening and nebulous background (Paper 1). In 
order to account for this variable background in the source 
detection, by following Paper 1, we estimated the complete- 
ness of two representative sub-regions i.e., regions A and C. 
90% completeness limits of these two regions have been ob- 
tained as 14.0, 13.5, 12.5 and 11.5 mag, respectively for 3.6, 
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4.5, 5.8 and 8.0 jim bands and are shown by using vertical 
dashed lines in Fig. \5\ 

The Sh2-252 complex lies at a distance of 2.4 kpc. The 
completeness limit of 14.0 mag at 3.6 /im corresponds to 
an approximate ste llar mass of 1.0 Mq for a YSO having 
an age of ~ 1 Myr (jSiess et al.l 120001 ). A more conservative 
estimate for the 90% completeness of the PMS membership 
of the region can be made from the histogram distribution 
of candidate YSOs shown by using the dashed lines in Fig. 
[S] which seems to be matching with the completeness es- 
timation by artificial star method. By taking variable ex- 
tinction and stellar crowding characteristics into account, in 
conclusion, we consider that our YSO census of the Sh2-252 
complex is complete down to 1-2 Mg. 



3.5 Characteristics of YSOs 



3.5.1 From SED an 

We constructed the spectral energy distribu tions (SEDs) us- 
ing the grid of models and fitting tools of iRobitaille et al.l 
l|2006l , 120071 ) for characterizing and understanding the na- 
ture of YSOs in the Sh2-252 compl ex. The models were 
computed by IRobitaille et al.l (|2006r ) by using a Monte- 
Carlo based 20,000 2- D rad iation transfer models from 
IWhitnev et"afl (|2003al lbl. [2004l ') and by adopting several com- 
binations of central star, disk, in-falling envelope, and bipo- 
lar cavity for a reasonably large parameter space. Each YSO 
model has SEDs for 10 viewing angles (inclinations), so the 
total YSO grid consists of 200,000 SEDs. The SED fitting 
tool provides the evolutionary stage and physical parame- 
ters such as disk mass, disk accretion rate and stellar tem- 
perature of YSOs and hence is an ideal tool to study the 
evolutionary status of YSOs in star forming regions. Of the 
total YSOs identified in Sections l3.3.1l and l3.3.2l 183 sources 
(51 Class I, 132 Class II) have counterparts in 24 /^m and 
we fit SEDs to these 183 YSOs. The counterparts of these 
YSOs are searched in the optical and NIR bands irrespective 
of their error criteria. Of these sources, 54 are found to have 
counterparts at least in two of the optical bands (UBVRI) 
along with the JHK, IRAC and 24 ^m data. 96 sources have 
JHK, IRAC and 24 /jm data and rest of the sources have 
IRAC and 24 /im data. Thus we have a minimum of 5 and 
maximum of 13 data points in the wavelength range from 
0.36 /urn to 24 jam to construct the SEDs. The SED fitting 
tool fits each of the models to the data allowing the distance 
and external foreground extinction as free parameters. We 
gave a distance range of 2.1 - 2.7 kpc (see Table 3 of Paper 
1) as the input. An extinction range has been supplied as 
the input for each source from the Ak value estimated from 
the (H — K) colour as discussed in Section T3. 2 1 and by allow- 
ing an uncertainty of ± 3 mag in the calculated extinction 
value. We further set photometric uncertainties of 10% for 
the optical and NIR data, 15% for the IRAC and 24 ^m 
data. These values are adopted instead of the formal errors 
in the catalog in order to fit without any possible bias caused 
by underestimating the flux uncertainties. We obtained the 
physical parameters of YSOs using the relative probability 
distribution for the stages of all the 'well-fit' models. The 
well-fit models of each source are defined by 



Xmin 



^ 2Ndata, 



where Xmin is the goodness-of-fit parameter for the 
best-fit model and Ndata is the number of input observa- 
tional data points. 

In Fig. |6j we show the sample SEDs of Class I, Class II 
and transition disk sources, where the solid black line rep- 
resents the best-fit and the grey lines are the subsequent 
well-fits. From the well-fit models for each source derived 
from the SED fitting tool, we calculated the x 2 weighted 
model parameters such as the stellar mass (M 4 ), temper- 
ature (T*), stellar age (£*), mass of the disk (Mdisk), disk 
accretion rate (Mdisk), envelope accretion rate (M env ) etc. 
The error in each parameter is calculated from the stan- 
dard deviation of all well-fit parameters. The parameters 
and the corresponding errors of all the YSOs are listed in 
Table [U The IDs 1 - 51 and 52 -183 represent Class I and 
Class II sources, respectively. For some of the sources, the 
errors associated with few parameters (Table [lj are quite 
large because we are dealing with a large number of param- 
eter space, with limited number of observational data points. 
Additional observational data points in longer wavelengths 
would help constrain these parameters more precisely. 

Table[T]shows that ~ 94% (172/183) of the sources have 
an age range between 0.1 to 5 Myr. This age distribution is 
in agreement with the age range discussed in Paper 1 using 
the optical colour magnitude diagram (CMD) of the candi- 
date YSOs. Similarly, masses of the YSOs are distributed 
in a range between 0.4 to 5.3 Mq with a majority of them 
having masses between 0.5 to 3.0 Mq. This is again in agree- 
ment with that of Paper 1. The SED model parameters of 
YSOs, such as the disk accretion rate and envelope accre- 
tion rate are considered to be functions of their evolutionary 
status. Stage 0/1, II and III YSOs have significant infalling 
envelopes, optically thick disks a nd optically thin disks, re- 
spectively (|Robitaille et al.ll2006P l. In Table [l] a majority of 
the Class II YSOs have the disk accretion rate of the order of 
~ 10 -7 - 10 -8 Mq yr _1 . Similarly, the envelope infall rates 
of most of the Class I YSOs are > 10 -6 Mg yr _1 , which are 
relatively higher than that of the Class II YSOs, confirming 
their evolutionary stages. In Fig. [7] we have shown the cumu- 
lative distribution of Class I and Class II YSOs as a function 
of their ages which manifests that Class I sources are rela- 
tively younger than Class II sources. We have performed 
the Kolmogorov-Smirnov (KS) test for the age distribution 
of the Class I and Class II sources. The test result shows 
that the chance of the two populations having been drawn 
from a same distribution is < 1%. It suggests that the age 
distribution of two populations are different. 



3.5.2 From optical colour-magnitude diagram 

The distribution of YSOs on the V/(V - I) CMD is an ideal 
tool to estimate their approximate ages and masses. The 
V/(V - I) CMD for the candidate Class I, Class II and Ha 
emission line sources are shown in the Fig. 8 of Paper 1, 
which suggests that majority of the YSOs in Sh2-252 are 
distributed between the PMS isochrones of age 0.1 and 5 
Myr and evolutionary tracks of masses ranging from 0.5 to 
3.0 Mq. This is in agreement with the SED based anal- 
ysis (Section I3.5.ip . The age of each YSO has been esti - 
mat ed from the CMD as discussed in lPandev et al.l (|2008f ) 
and lChauhan et al.l (|2009l ). Briefly, their ages have been de- 
rived by comparing their locations on the CMD with PMS 
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Figure 5. Apparent magnitude histograms of detection in all the four IRAC bands having errors < 0.2 mag are shown with solid 
histograms and dashed histograms are for the YSO detection. The vertical dashed lines indicate 90% completeness limit estimated from 
the artificial star method. 
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Figure 6. Sample SE Ds for a Class I (left), Class II (middle) and transition disk (right) sources using the SED fitting tools of 
iRobitaille et al.l |2007]). The solid black line shows the best-fit and the grey lines show the subsequent well- fits. The dashed line shows 
the stellar photosphere corresponding to the central source of the best fitting model. The filled circles denote the input flux values. 
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Figure 7. Cumulative distribution of Class I (solid) and Class II 
(dotted) YSOs as a function of age. 
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Figure 8. Comparison of age distribution from CMD analysis 
with that from SED fitting. Squares and tringles represent the 
Class I and Class II YSOs, respectively 



isochrones of various ages (|Siess et alJl200Cl ) , after correcting 
for the distance and reddening of the region. Apart from the 
Class I, Class II and Ha sources, all other sources having age 
< 5 Myr in the V/(V-I) CMD and with {H-K)> 0.4 mag 
were suggested as the candidate PMS members of the com- 
plex in Paper 1. We also estimated the ages of these addi- 



tional candidate PMS sources as discussed above. There are 
80 YSOs which have age estimation both from the CMD and 
SED (Section I3,5.1[l analyses. In Fig. [8] we have shown the 
comparison between the age distribution from CMD analysis 
with that from SED fitting. Squares and triangles represent 
the Class I and Class II YSOs, respectively. Though the age 
estimation from both the methods are quite scattered, in 
general, the distribution can be considered correlated with 
each other. We have also performed the KS test between 
the age distribution of YSOs from SED and CMD analyses 
and found that the ages estimated by these two methods are 
correlated by a factor of ~ 81%. Hence, we confirm that our 
age estimation based on both methods is in agreement. 



3.6 The ages of sub- regions in Sh2-252 

We attempted to constrain the evolutionary status of the 
sub-regions of Sh2-252 using the SED fitting and CMD anal- 
yses. Here we would like to caution that although the age 
estimation using the SED fittings may not be robust, the 
results agree reasonably well with the ages estimated by us 
using V/(V - I) CMD for a subset of YSOs (Section [3X2). 
Hence we accept the ages determined by models with rea- 
sonable confidence. We estimated the ages of the sub-regions 
A, C, E, NGC 2175s and Teu 136 of Sh2-252 by using the 
ages of the candidate YSOs lying within their boundaries. 
The boundaries of these individual regions have been taken 
from the stellar surface density map. i.e., 3'. 5, 3'. 5, 3'.0, 3'.0 
and 3'.0, respectively for regions A, C, E, NGC 2175s and 
Teu 136 (Paper 1). The histograms of age distribution esti- 
mated from the SED and CMD analyses of each sub-region 
are shown in Fig. [5] The median ages of YSOs within A, C, 
E, NGC 2175s and Teu 136 have been obtained as ~ 0.5, 
1.5, 1.5, 2.5 and 2.5 Myr, respectively. As evident in the 
K-band surface stellar density map (Paper 1), no cluster- 
ing is apparently seen towards the region B. From the SED 
analysis of 3 candidate YSOs lying within 2'. 5 radius of re- 
gion B, a median age of 1-2 Myr can be assigned for this 
region. However, the age of Sh2-252 B should be considered 
as an approximate value since the statistics is very poor. 
The above age estimates based on the SED/CMD analyses 
show that the region A is the youngest among all the sub- 
regions, whereas regions B, C and E are of similar age and 
NGC 2175s and Teu 136 are found to be slightly evolved 
than other sub-regions. 



8 3.7 Spatial distribution of YSOs in Sh2-252 

Spatial distribution of YSOs in a star-forming complex helps 
to trace its star formation history. In Fig. 1101 we show the 
spatial distribution of all the identified candidate YSOs in 
Sh2-252 (i.e., red squares: Class I; green triangles: Class 
II; magenta circles: transition disk sources) overlaid on the 
IRAC 4.5 /im mosaic image. The 61 Ha emission line sources 
identified from slitless spectroscopy survey (Paper 1) are 
shown by using blue crosses. The important sub-regions of 
Sh2-252 are also marked in the figure. This figure indicates 
some interesting spatial sequence of YSOs within Sh2-252. 
A comparison of this figure with Fig. [5] clearly indicates 
that a majority of the candidate YSOs identified in Sh2- 
252 are preferentially concentrated around the CH II regions 
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Figure 9. Histograms of age distribution estimated from the SED and CMD analyses for candidate YSOs within the five sub-regions of 
Sh2-252 complex with median age of 0.5, 1.5, 1.5, 2.5 and 2.5 Myr, respectively for the regions A, C, E, NGC 2175s and Teu 136. 



A, C and E, respectively. Also, the bright rim feature dis- 
cussed in Section [T] (i.e., region F), which makes up the 
western border of the complex contains several number of 
YSOs. There are not many YSOs detected towards the re- 
gion B. The WCF is seen to have more number of YSOs 
when compared to the ECF. Overall, the distribution of 
YSOs in the west of Sh2-252 shows a nice correlation with 
the semi-circular shaped distribution of molecular gas and 
the PAH emission/H II region cavity boundaries (see Section 
13. l| l. The strong positional coincidence between the YSOs 
and the molecular cloud towards the WCF in Sh2-252 shows 
the enhanced star formation activity towards this region 
as observed in other s tar forming regions (e.g. I Evans et al.l 
12009 . iFang et al.ll2009l ). There are a few number of sources 
scattered towards the east of Sh2-252. There are two groups 
of YSOs aligned in a filamentary manner associated with 
the cloud remnants seen towards the east of CH II region E 
(see Figs, [l] and [3}. However, the cluster NGC 2175s does 
not seem to have significant number of YSOs. There seems 
to be a small clustering of Class II and transition disk YSOs 
towards Teu 136. It is also interesting to note that there is 
an enhanced concentration of Class I YSOs towards the east 
of the sub-region A when compared to other sub-regions of 
Sh2-252. This is an indication that the region A might be at 
a younger evolutionary stage and we will be iscussing more 
on this in the ensuing sections. 



The number of YSOs in each class is use ful for estimat- 
ing the relative age of a star forming region (|Hatchell et al.l 
120071 ; iGutermuth et alj|2009f ). Since Class II sources are as- 
sumed to be more evolved than Class I and hence the Class 
I/II ratio is a function of the evolutionary status of a sys- 
tem (|Beerer et al.l 120101 1. We used this ratio as a proxy for 
the relative ages of the sub-regions in Sh2-252. The spatial 
distribution of the candidate YSOs given in Fig. [TU] shows 
that most of them are distributed within or in the close 
proximity of the sub-regions A, C, E and Teu 136. Hence we 
estimated the ratio of Class I to Class II YSOs for these sub- 
regions within their radius (Paper 1). The Class I/II ratios 
have been obtained as 0.54, 0.24, 0.36 and 0.13, respectively 
for regions A, C, E and Teu 136. This fraction is of almost 
similar order with the median ratio of C lass I to Class II 
sourc es (0.27) obtained in nearby clusters (|Gutermuth et al.l 
120091 ). The various sub-regions within the S254-258 complex, 
a nearby star forming region to Sh2-252 in the Gem OBI as- 
sociatio n, also have the Class I/II ratios vary between 0.32 
to 0.78 (|Chavarrfa et all 120081 ) . Thus the Class I/II ratios 
of the sub-regions of Sh2-252 can be considered comparable 
with the other star forming regions. In this scheme, the high 
concentration of Class I YSOs towards the region A suggests 
that the youngest population is located towards the western 
boundary of the complex. 
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Figure 10. Spatial distribution of candidate YSOs in Sh2-252 over plotted on the 4.5 fira mosaic image. The candidate Class I, Class II and transition disk sources are shown by using 
red squares, green triangles and magenta circles, respectively and the blue crosses are the Ha emission line sources identified from the slitless spectroscopy survey. The main sub-regions 
of Sh2-252 are marked in the figure and the star symbol represents the location of HD 42088. RA and Dec are in J2000 epoch (see the electronic version for the colour image). 
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4 DISCUSSION 

4.1 Probable evolutionary status of the 
sub-regions in Sh2-252 

Since most of the giant molecular clouds (GMCs) in the 
Galaxy are hierarchical in structure, the evolutionary sta- 
tus of the identified YSOs helps to get a deep insi ght about 
the ov erall star formation activity of the complex. iTei et al.l 
([20061 ) studied the sub-region A using NIR, sub-mm and 
radio observations. Based on their K-hemd surface density 
map, which is deeper and of better resolution t han ours, 
they f ound sub-structures within this region. In ITei et al.l 
1)20061 ) . the prominent peak of the stellar density coincides 
spatially with the central ionizing source of A, which is an 
09.5V type star (Paper 1). However, a secondary stellar 
density peak, which lies ~ 50" east to the main peak has 
also been found and its spatial locat ion coincides wi th the 
sub-mm peak detected at 850 fj,m (|Tei et all [20061 ) . The 
authors have attributed this as a possible cold core at an 
early evolutionary stage than A. Similarly, the most intense 
CO peak of the WCF is also located east to the region A 



with water and methanol maser emissions ( Lada fc Wooden! 
1 19791 ; iKoempe et"ai]|l989l ; ISzvmczak et al]|200Ch within its 



proximity, suggesti ng the recent star formation activity to- 
wards the east of A llLada fc Woodenlll979l ; lLada et al]|l98ll ; 
IKoempe et al.lll989T ). The water and methanol maser emis- 
sions associated with 850 /im and CO peak towards the east 
of region A suggests that this region is possibly at an early 
evolutionary stage than A. 

The high concentration of Class I YSOs around the re- 
gion A (Fig. I10|) and the large value of Class I/II ratio (see 
Section I3.7f) show that region A is at a young evolution- 
ary stage compared to the other sub-regions of the complex. 
The average age estimation (~ 0.5 Myr) based on the SED 
and CMD analyses also shows that the region A is indeed 
very young. Based on the class I/II ratio and age histogram 
analysis, the regions B, C and E are found to be at a similar 
evolutionary stage and are slightly evolved (median age ~ 
1-2 Myr) than region A. The presence of a very few number 
of YSOs towards the small cluster NGC 2175s (see Fig. 1X0)1 
suggests that either the cluster was formed in an early epoch 
of star formation in Sh2-252 or it i s a case of rapid circum- 
stellar disk evolution of YSOs fe.g jBallv et alj|l998h due to 
the presence of four early type massive sources at its center 
(see Table 4 of Paper 1). However, the median age estima- 
tion based on the CMD analysis gives an approximate age of 
~ 2-3 Myr for this cluster. The small cluster Teu 136 lying 
towards the east of Sh2-252 seems to have more number of 
Class II YSOs than Class I. The low value of the class ratio 
as well as the SED/CMD based age analyses (median age 
~ 2-3 Myr) show that this cluster is also at a slightly more 
evolved stage when compared to the sub-regions A, B, C, 
and E of the Sh2-252 complex. 



4.2 Star formation activity towards WCF-The 
general picture 

Sh2-252 is a large H n region in the Gem OBI association, 
excited by the 06.5V star HD 42088 of age 2-4 Myr. It 
contains 12 B-type stars, some of them are associated with 
CH n regions. Since CH II regions are at a early stage of 



massive star formation, their impact is limited to their lo- 
cal environment. Thus, 06.5V star is the most dominant 
source of energy respo nsible for ionization in the region. 
lLada fc Wooden! |l979j) showed that the cloud complex is 
separated mainly into two fragments ECF and WCF, which 
themselves are embedded in a more diffuse and extended (~ 
90 pc) molecular cloud complex. Most of the mass of WCF 
(~ 2.5 x 10 4 M©) is concentrated in a narrow ridge, which 
borders the western IF of Sh2-252. This dense, semi-circular 
shaped shell observed at the eastern boundary of WCF, ap- 
pears to have been swept up by the pressur e from the IF gen- 
erate d by the exciting source of Sh2-252 (JLada fc Wooden! 
ll979T ). The large size and massive content of WCF suggest 
that the entire cloud is unlikely to be swept up. We roughly 
compared the required kinetic energy and the energy pro- 
duced by the massive 06.5V star to sweep up the entire 
shell. The combined energy (radiation plus wind) liberated 
by the 06.5V stars is found to be less than the required en- 
ergy to drive the entir e cloud (WCF+ECF ~ 2.7xl0 4 M ; 
lLada fc Wooden! 1 19791 ). Thus it is suggested that at least 
part of the shell in WCF is probably primordial in nature 
and the entire shell is not swept up by the expansion of the 
H II region. However, at the same time it is obvious from 
Fig. [5]that the H II region is in interaction with WCF, with 
bright semi-circular PDR and compressed neutral material 
behind it and protruding structures such as 'finger tips' and 
a bright rim pointing towards the ionizing source. More- 
over, we observed recent star formation activity towards this 
region. The morphology resembles that star formation has 
occurred in thin shells around the H II region and star for- 
mation seems to be more active in this direction compared 
to the whole cloud complex. 

Fig. [TT] shows the colour composite image of the west- 
ern half of Sh2-252 made from the 8.0 and 24 ^m bands. 
The contours of 1.1 mm dust contin uum emission (taken 
from Bolocam Galactic Plane Survey; lAguirre et al.l 1201 1! ) 
and low resolution radio continuum emission at 1280 MHz 
(Samal et al. in preparation) are over plotted in Fig. 1111 
On a large scale we see a ring of molecular material at 1.1 
mm aligned parallel to the outer extent of the IF of Sh2- 
252 seen in optical (see Fig. [1} , and located just behind the 
bright PAH emission seen in 8 fim and also corresponds to 
the dark region at 24 /an. The molecular material in the 
ring seems to be fragmented at various locations. This shell 
of molecular material located behind the IF provides evi- 
dence for the most compressed zone of the WCF fragment. 
We also observed many Class II and Class I YSOs in this di- 
rection (Fig. I10p . with a majority of them projected behind 
the PDR, in the vicinity of 1.1 mm emission. Apart from the 
YSOs, three CH II regions (A, B and C) are observed in the 
close vicinity of the 8 /im ring, with region A being located 
slightly behind the PAH ring and separated from the ring by 
a dark lane, consisting of highly dense molecular material. 
The region A is part of the massive (M ~ 2.5xl0 4 M©), 
large (~ 16 pc) WCF cloud (see Fig.[T}. Hence the influence 
of 06.5V star to region A, if at all exists, must be weak. 

In general, the configuration (at least in the western di- 
rection) suggests that the H II region has collected some of 
the matter during its expansion at its periphery, like the col- 
lect fc collapse process of star formation such as in Sh2-212 
(|Deharveng et alj|2004 1 Jose et al.ll201ll ). Collect & collapse 
process works well in homogeneous and uniform medium, 



Young stellar population in Sh2-252 15 




Figure 11. Colour composite image of the western half of Sh2-252 made from the 8.0 and 24 fim bands along with the 1.1 mm dust 
continuum emission map (white; contour levels at 0.06, 0.21, 0.37, 0.52, 0.68, 0.83, 0.99, 1.14 and 1.30 Jy/beam) and the low resolution 
radio continuum map at 1280 MHz (green; contour levels at 0.25, 0.44, 0.63, 0.82, 1.02 mjy/beam). The sub-regions are marked in the 
figure. North is up and east is to the left. 



which is definitely not the case here, with openings at various 
directions and without a clear circular shell structure. How- 
ever it has been found that collect & collapse process can sti ll 
be applicable in turbulent medium (jDeharveng et al.ll2008h . 
From the extinction map (Fig. [3]) we found that the visual 
extinction in the direction of the PAH ring is, Ay > 5.5 
mag, which corresponds to a column density of hydrogen, 
iV(H2) in the region > 1.03xl0 22 atoms cm" 2 , using the re- 
lation A v = 5.34 xlCT 22 iV(H 2 ), derived from iBohlin et all 
1 19781 ') for A v = 3.1xE(B-V). The observed column den- 
sity towards the molecular ring is thus comparable to the 
column density required (i.e., ~ 6xl0 21 atoms cm -2 ) in the 
co llected layer for the fragm entation to happen, as predicted 
bv lWhitworth et al.l 1)19941 ). Assuming the observed star for- 
mation near the vicinity of the IF might have happened 
due to collect & collapse process, we compared the observed 
pr operties with the predica tion by the numerical simulation 
of IWhitworth et al.l (|l994h . We ad opted 1.7xl0 49 as ion- 
izing photons s _1 (for 06.5V star; IVacca et al.lll996T ). the 
isothermal sound speed as 0.2-0.3 km s~ and the density of 
the neutral material into which t he H n region evolved as 
400 cm" 3 (|Lada fc Wooded 1 1979h . With these values, the 
predicted fragmentation time and radius at which fragmen- 
tation would take place are found to be 2.3-2.9 Myr and 
10.3-11.6 pc, respectively. Though there are large uncertain- 
ties involved in these calculations, considering the age of the 
ionizing star is of ~ 2-4 Myr (Paper 1) and ages of the Class 
I/II YSOs projected on the shell ~ 0.5 Myr, as well as the 
fact that the farthest separation of the shell from the ioniz- 



ing star is ~ 10 pc, the star formation activity could in fact 
might have taken place in the collected shell. 

The possible formation of regions B, C and E is in 
doubt by the same scenario, owing to their large age (~ 1-2 
Myr; Section [4T}. These objects might have formed before 
the ionization-shock front passed over them, although with 
the present data it is hard to prove whether the passage of 
ionization-shock front makes the age of these sources look 
older. Ignoring this fact, we suspect that they are possibly 
not triggered due to the expanding H n region, and more 
likely to be primordial in nature, resulted from the forma- 
tion and evolut ion of the molecu l ar clo ud itself, which is in 
agreement with lLada fc Wooden! [1 19791 ). 

The radio source F is an extended one, with no inter- 
nal heating source. The radio emission is located ahead of 
the 8 /im emission (see Fig. [T5] for a close up view of region 
F) and corresponds to the glowing gas in the optical image 
(see Fig. 2 of Paper 1), thus, more likely due to the photo- 
evaporating gas from the surface of the cloud exposed to 
radiation. This phenomenon usually occurs when a massive 
star drives an ionization-shock front into the ISM, the ex- 
ternal layers of the cloud exposed to radiation starts photo- 
evaporating, and can be seen as glowing layer of gas. In Fig. 
[T3l the HST-WFPC2 archive image at 0.65 ^m (Program 
No: 9091; PI: J. Hester) of the extreme edge of the molecular 
ridge at region F (i.e., the brightest portion of the bright- 
rim feature in Fig. 1 of Paper 1) is shown. Apart from the 
photo-evaporating gas, this better resolution image displays 
number of 'finger-tips and pillars' like fragments which are 
seen silhouetted against the bright background pointing to- 
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Figure 12. Close-up colour composite image of region F made 
from the 8.0 and 24 (im bands along with the 1.1 mm dust con- 
tinuum emission map (white) and the low resolution continuum 
map at 1280 MHz (green). North is up and east is to the left. 



wards the direction of HD 42088. Three YSOs are observed 
at the tip of region F (see Fig. I10[) and they are low mass 
(< 2M0) Class I/II sources, thus possibly formed from low 
mass clouds. Assuming that the impact of ionizing radiation 
to such cloud structures initiates new star formation similar 
to those in BRCs and small globules, we compared the ob- 
served incident ionizing flux to the required range of ionizing 
fluxes tha t can trigger s t ar for mation in fiducial cloud as pre- 
dicted by Bisbas et al.l (|2011T ) using numerical simulations. 
iBisbas et al.l ( 20 111 ) numerically estimated the required ion- 

3xlO n cm -2 s _1 . They 



izing flux range as 10 9 cm -2 s~ x 
predicted that if the ionizing flux is very low (<10 9 cm -2 
s _1 ) then the cloud will be slowly eroded and if the ionizing 
flux is high (>10 11 cm~ 2 s _1 ) then the cloud will be rapidly 
dispersed, with no star formation. We calculated the ion- 
izing photon flux impinging onto these structures from the 
radio emission from the surface layer of the cloud following 
iLefloch et al.l (J1997T ) assuming the electron temperature of 
the photo-evaporating gas as 10 4 K, and it is estimated to 



be 3xlO a cm" 



The observed estimation is well within 



the predicted range to initiate star formation. 

Besides the low-mass YSOs at the tip of the pillars, we 
also observed a clear clustering or concentration of YSOs in 
the dense region of the cloud located behind region F (see 
Figs. Hand [TUJ). The age of these YSOs is of the order of ~ 
1 Myr, thus indicates the youth of the region. These YSOs 
seem to have resulted from the gravitational fragmentation 
of the collected material rather than the implosion of pre- 
existing clumps due to external pressure. The YSOs show 
a chain-like elongated distribution parallel to the IF and 
appear to lie in the collected shell. Whereas in the case of 
radiation driven implosion (RDI), the YSOs are expected 
to lie perpendicular to the direction of progressive IF, with 
younger YSOs being located away from the IF. The incident 
photon flux is also on the lower-side of the flux needed to 
initiate RDI, thus possibly is not enough for a massive cloud 



to form a cluster of stars. Hence, we suspect that the star 
formation is more likely due to gravitational instability in 
the collected material, such as seen in the massive condensa- 
tion of RC W120, where a chain o f YSOs is observed parallel 
to the IF (|Deharveng et alj|2009h . 

In Fig. \W\ we find the densest concentration of YSOs 
in the close vicinity of the molecular clump lo cated be- 
tween the CH 11 regions A and B (see Fig. 15 of lTei et al.l 
I2006J ). with many YSOs of Class I nature (~ 10 5 yr). This 
cloud clump is associated with water and methanol masers 
and do es not show any radio emission (down to 0.4 mjy). 
Hence, iTei et al.l (|2006l ) suggested that it could be a site 
of high-mass star-forming protocluster in a very early evo- 
lutionary stage (Section 14.10 . We also observed a cluster 
of YSOs, mostly Class I (Fig. I10p . thus reflecting it is in- 
deed a site of cluster formation sandwiched between the two 
relatively evolved CH 11 regions A a nd B. Using H2 obser- 
vations at 2.12 urn. ITei et alj (|2006l ) traced a cometary-arc 
shaped shocked molecular H2 emission bounding the ionized 
gas (i.e., ionization bounded side of the CH 11 region A) in 
the direction of the molecular clump, possibly reflecting the 
interaction of the H 11 region with the adjacent molecular 
cloud. Similarly, we noticed sharp radio contours for region 
B, in the direction of the molecular clump (Fig. Ill|) . The 
radio size of region B is larger than regions A, C and E, and 
is ionized by a similar kind of star (B0-B0. 5 type; Paper 1). 
Thus, if they are evolved in a cloud clump of similar den- 
sity, the larger size of region B possibly reflects that it is 
relatively evolved in nature. The region between A and B 
indeed reflects a site of recent star formation, where much 
younger sources are found in the matter in between the two 
relatively evolved CH 11 regions A and B. A similar kind of 
star formation has been observed at the interfaces of two ex- 



panding super giant shells LMC4 and LMC5 (I Cohen et al 
200S ) and two classical H 11 regions S255-S257 (|Oiha et al 



20111 ). where it is believed that collision of dense swept up 
neutral materi als associated with indiv idual shells lead to 
star formation (|Yamaguchi et al.ll2001al lbl). Possibly, we are 
witnessing one such region, where star formation has been 
enhanced due to a similar kind of process. 

The overall picture emerges from our analyses is that 
the star formation activity observed towards WCF is more 
likely multi-fold and the sources are found to be at differ- 
ent evolutionary stages . Though the large scale CO map 
IjLada fc Wooden! 1 19791 ) provides much information about 
the global structure of the complex, but the region certainly 
deserves attention for high resolution molecular line obser- 
vations to explore the detailed properties of the sub-regions. 



5 SUMMARY 

The optically bright, extended H 11 region Sh2-252 is a part 
of the Gemini OBI association. This region is mainly com- 
posed of two small clusters NGC 2175s and Teu 136 and 
four CH 11 regions namely A, B, C and E. In this paper, an 
extensive survey of the star forming complex Sh2-252 has 
been undertaken with an aim to explore its hidden young 
stellar population, their characteristics, spatial distribution, 
morphology of the region and finally to understand the star 
formation scenario of the complex for the first time. Spitzer- 
IRAC, MIPS photometry (3.6 - 24 fj,mj are combined with 
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Figure 13. HST-WFPC2 archive image at 0.65 /im of the sharp edge of the molecular ridge in region F. The direction of the ionizing 
source HD 42088 is marked in the figure. 



2MASS-NIR and optical data sets to identify and classify 
the YSOs by their IR excess emission from their circumstel- 
lar material. Using the IR C-C criteria, we have identified 
577 YSOs in the complex, of which, 163 are consistent with 
Class I, 400 are consistent with Class II and 14 are con- 
sistent with transition disk YSOs, suggesting a moderately 
rich number of YSOs in this H II region. From the CMD 
and SED based analyses, majority of the YSOs are found to 
have an age distribution between 0.1-5 Myr and mass in 
the range of 0.3 - 3.0 M . 

Spatial distribution of the candidate YSOs shows that 
they are mostly clustered around the sub- regions of the com- 
plex such as in A, C, E, F and Teu 136. However, majority 
of the candidate YSOs are distributed in the western part 
of the complex when compared to the east, suggesting en- 
hanced star formation activity towards its west. Using the 
SED and CMD based age analyses, we derived probable evo- 
lutionary status of the sub-regions of Sh2-252. Our analyses 
suggest that the region A is the youngest (~ 0.5 Myr), the 
regions B, C and E are of similar evolutionary stage (~ 1-2 
Myr) and the small clusters NGC 2175s and Teu 136, lo- 
cated towards the east of Sh2-252 are slightly evolved (~ 
2-3 Myr). Morphology of the region in the Spitzer colour 
images as well as in the 1.1 mm map shows an almost semi- 
circular ring-like shape towards the western half of the com- 
plex. Indeed, we find a molecular shell composed of several 
clumps distributed around the main ionizing source (i.e., HD 
42088), suggesting that the expansion of the H II region is 
collecting the surrounding material, which gives rise to the 
semi-circular ring shape. We find several candidate YSOs 
distributed over the semi-circular molecular shell which is 



an evidence for the star formation activity within the shell. 
Finally, by comparing the age of the ionizing source, frag- 
mentation time of the collected molecular shell and age of 
the YSOs, we suggest collect and collapse scenario as one of 
the the possible mechanisms responsible for the star forma- 
tion within the shell. We observed the densest concentration 
of YSOs, (mostly Class I, ~ 0.5 Myr) at the western out- 
skirts of the complex, within a molecular clump located be- 
tween the CH n regions A and B. The correlation between 
the molecular clump at this location with large number of 
Class I YSOs and the associated water and methanol masers 
suggest that it is indeed a site of cluster formation in a very 
early evolutionary stage, sandwiched between the two rela- 
tively evolved CH II regions A and B. We conclude that the 
region is undergoing a complex star formation activity and 
there is a strong interplay between the radiation from the 
expanding H II region to the surrounding molecular mate- 
rial. The region certainly deserves attention for high resolu- 
tion molecular line observations to further explore its hidden 
structure. 
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Table 1. Physical parameters of YSOs from SED fits. The complete table is available in the electronic version. 



ID 


"(2000) 


•5(2000) 


M, 


T, 


t» 


Mdisk 


Mdisk 


Me„ 




(deg) 


(deg) 


(*%) 


(10 4 K) 


(10 6 yr) 


(M & ) 


(10- 7 A/ o /yr) 


(10- 6 A/ Q /yr) 


1 


92.242438 


20.391663 


2.1 ± 0.5 


0.77 ± 0.29 


4.42 ± 3.28 


0.026 ± 0.028 


3.542 ± 2.896 


0.002 ± 0.003 


2 


92.066338 


20.544426 


1.4 ± 1.0 


0.43 ± 0.12 


0.06 ± 0.23 


0.019 ± 0.040 


5.907 ± 5.553 


126.0 ± 116.0 


3 


92.196279 


20.630317 


4.5 ± 0.3 


0.75 ± 0.10 


0.74 ± 0.23 


0.001 ± 0.001 


0.038 ± 0.037 


0.400 ± 0.562 


4 


92.464273 


20.476278 


3.1 ± 0.5 


1.17 ± 0.15 


4.85 ± 1.77 


0.003 ± 0.007 


0.028 ± 0.028 


0.000 ± 0.000 


5 


92.068404 


20.607574 


3.7 ± 0.6 


1.34 ± 0.14 


4.00 ± 2.17 


0.011 ± 0.022 


1.256 ± 1.186 


1.127 ± 17.51 



